The quasiclassical limit of the symmetry 
constraint of the KP hierarchy and the 
dispersionless KP hierarchy with 
self-consistent sources 

Ting Xiao Yunbo Zengf 

Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 

f Email: yzeng@math. tsinghua. edu. cn 

Abstract 

For the first time we show that the quasiclassical limit of the symmetry constraint 
of the KP hierarchy leads to the generalized Zakharov reduction of the dispersionless 
KP (dKP) hierarchy which has been proved to be result of symmetry constraint of 
the dKP hierarchy recently. By either regarding the constrained dKP hierarchy as 
its stationary case or taking the dispersionless limit of the KP hierarchy with self- 
consistent sources directly, we construct a new integrable dispersionless hierarchy, i.e., 
the dKP hierarchy with self-consistent sources and find its associated conservation 
equations (or equations of Hamilton- Jacobi type). Some solutions of the dKP equation 
with self-consistent sources are also obtained by hodograph transformations. 

Keywords: quasiclassical limit; constrained KP hierarchy; (dispersionless) KP hierar- 
chy with self-consistent sources; conservation equation (equation of Hamilton- Jacobi type); 
Zakharov reduction; hodograph transformation 

1 Introduction 

The dispersionless integrable hierarchies provide us an interesting type of nonlinear inte- 
grable models which have important applications from complex analysis to topological field 
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theory (see [1-12] and the references therein). In [1, 2, 5, 6, 7], a standard procedure of 
dispersionless limit of integrable dispersionfull hierarchies is proposed. In this procedure, 
dispersionless hierarchies arise as the quasiclassical limit of the original dispersionfull Lax 
equations performed by replacing operators by phase space functions, commutators by Pois- 
son brackets and the role of Lax pair equations by conservation equations (or equations 
of Hamilton- Jacobi type ). Ad scheme of dispersionless hierarchies has been proposed by 
Konopelchenko et al (see [11, 12] and the references therein). Recently, from this point 
of view, some important reductions of dispersionless hierarchies are shown to be nothing 
but symmetry constraints [12]. Also several methods for solving dispersionless hierarchies 
have been formulated such as twistorial method [7, 13] and hodograph transformation [5, 6]. 
In [5, 6], from the conservation equations of the dispersionless KP equation, Kodama and 
Gibbons found exact solutions of it and its reductions by hodograph transformations and 
obtained general hodograph equations for hydrodynamical type equations. 

The soliton equations with self-consistent sources (SESCS) are another type of impor- 
tant integrable models in many fields of physics, such as hydrodynamics, state physics, 
plasma physics, etc (see [14-18]). For example, the KP equation with self-consistent sources 
(KPESCS) describes the interaction of a long wave with a short-wave packet propagating on 
the x,y plane at an angle to each other (see [14] and the references therein). In general, the 
constrained integrable hierarchy may be viewed as stationary system of the corresponding 
hierarchy with self-consistent sources and the Lax representation for the latter can be ob- 
tained naturally from that for the former [16-18]. In this sense, the soliton hierarchies with 
self-consistent sources may be viewed as integrable generalizations of the original soliton 
hierarchies. 

In contrast with the dispersionfull integrable hierarchies with self-consistent sources, the 
dispersionless integrable hierarchies with self-consistent sources have not been studied be- 
fore. In this paper, we first investigate the quasiclassical limit of the symmetry constraint 
of the KP hierarchy which leads to the generalized Zakharov reduction of the dKP hierar- 
chy. The latter has recently been proved to be symmetry constraint of the dKP hierarchy 
using the 8 method [12]. So, we find the relation that the symmetry constraint for the Lax 
function of the dispersionless hierarchy can be obtained by the quasiclassical limit of the 
symmetry constraint for the Lax operator of the corresponding dispersionfull hierarchy. By 
either regarding the constrained dKP hierarchy as its stationary case or taking the disper- 
sionless limit of the KP hierarchy with self-consistent sources directly, we construct a new 
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dispersionless hierarchy, i.e., the dKP hierarchy with self-consistent sources. This hierarchy 
is also integrable since we can find its associated conservation equations (or equations of 
Hamilton- Jacobi type). So in this sense, the dKP hierarchy with self-consistent sources may 
be viewed as an integrable generalization of the dKP hierarchy. Compared with the case 
of a nonlocal term appearing in the t-part of the Lax pair for the KP equation with self- 
consistent sources [17], the t-part of the conservation equations for the dKP equation with 
self-consistent sources possesses rational terms with poles. From the obtained conservation 
equations, we can solve the dKP hierarchy with self-consistent sources by hodograph trans- 
formations. 

The paper will be organized as follows. In section 2, we remind some definitions and 
results about the KP hierarchy with self-consistent sources. In section 3, we take the qua- 
siclassical limit to the symmetry constraint of the KP hierarchy to obtain the symmetry 
constraint of the dKP hierarchy. The dKP hierarchy with self-consistent sources and its 
associated conservation equations are obtained. Section 4 is devoted to solving the dKP 
hierarchy with self-consistent sources by hodograph transformations and some solutions for 
the dKP equation with self-consistent sources are presented. 

2 The KP hierarchy with self-consistent sources 

We first review some definitions and results about the KP hierarchy with self-consistent 
sources in the framework of Sato theory [17-20]. Given a pseudo-differential operator (PDO) 
of the form 



where d = d x , dd 1 = d 1 d — 1, t — (ti — x, t 2: • • • ), the KP hierarchy is defined as 



where B n = (L n )> is the differential part of L n . The Lax euqation (2.2) is equivalent to the 
existence of the Baker function -0 such that 



L = d + u 1 (t)d~ 1 +u 2 {t)d~ 2 + 



(2.1) 



d tn L = [B n ,L], 



(2.2) 



Lip = \tp, 



(2.3a) 



d tn i) = B n ip, 



(2.3b) 



and ip also satisfies 



d x ij} = Mip, 



(2.4) 
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where M is the Orlov operator of the KP hierarchy. The adjoint Baker function ifj* satisfies 

L*i)* = Xif;*, (2.5a) 

d n r = -BW, (2.5b) 
d x tp* = -M*tj)*. (2.5c) 
Making a constraint of the PDO L (2.1) as [21] 
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L n = (L")> + Y^q t (t)d- 1 r l (t), n G N (2.6) 
l 

where and rj(i) satisfying 

%t m = B m q h r ittm = -B* m n, B m = [(L n )f ]> , z = 1, ...,7V, Vm e N, (2.7) 

we will get the n-constrained KP hierarchy as 

(L B ) tm = [S mj L B ] J (2.8a) 

%t m = B m qi, (2.8b) 

r htm = -B* mri , i = l,...,N. (2.8c) 

If we add the term (B m ) tn to the right side of equation (2.8a), the KP hierarchy with self- 
consistent sources will be obtained as [17] 

(B m ) tn -(L n ) tm + [B m ,L n ]=0, (2.9a) 

qi,t m = B m qi, (2.9b) 

^ = -fl* m ri, i = l,...,N. (2.9c) 

As many cases in (1+1) -dimension (see [16] and the references therein), the n-constrained KP 
hierarchy may be considered as the stationary one of the KP hierarchy with self-consistent 
sources if "t n " is viewed as the evolution variable and the Lax representation for the KP 
hierarchy with self-consistent sources can be obtained naturally from that for the constrained 
KP hierarchy [17]. The most important equation in the hierarchy is the KP equation with 
self-consistent sources obtained when n = 3, m = 2 in (2.9), [14] 

1 N 3 

[u 1;t - 3ui«i )X - -u l<xxx + ^2(qiTi)x]x ~ j u i,yy = 0, (2.10a) 
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q%, y = %xx + 2mqi, (2.10b) 

n, y = -r itXX - 2uin, i=l,...,N (2.10c) 

where t — t s ,y — t 2 . With (2.10b) and (2.10c), (2.10a) will be obtained by the compatibility 
of the following linear equations (Lax pair) [17] 

Vv = i>xx + 2wiV>, (2.11a) 

3 _ N 

ipt = i^xxx + 3wi^x + -K,* + (d 1 u hy ))il> + ^2q i d~ 1 (r i 'if)). (2.11b) 

i=i 

3 Dispersionless limit 

Following the standard procedure of dispersionless limit introduced in [1, 2, 5, 6, 7], we 
will get the dispersionless counterpart of (2.9) which can be regarded as the dispersionless 
KP hierarchy with self-consistent sources. Simply taking T n = et n and thinking of u n (j) = 
U n (T) + O(e) as e -> 0, L in (2.1) changes into 

oo „ oo 

L e = ed + J2 u i(-)(^y i = ^ + ^(^(T) + 0(e))(ea)- i , d = d x , X = ex. (3.1) 
i=i e i=i 

The constraint (2.6) now changes into 

N 

L: = B tn + Y, qt {-){edr l r t {-), 5 en = (L?)> , (3.2) 



e e 
i=i 

where ^(7) and ri(j) satisfy 

c[ft(|)]T m = B emft (^), c[ri(^)]T m = -S^l), S em = [(^)-]>o, < = 1,...,JV, (3.3) 
and the counterpart of (2.8) is 

e(^) Tm = [5 €m ,^], (3.4a) 
tk(j)}T m = £«„fc(|), (3.4b) 

c[ri(|)]r m = -iO-i(f). < = 1,...,JV. (3.4c) 

It is proved in [7] that 

00 

C = a e (L e )=p + Y,U i (T)p- i . (3.5) 
i=i 
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is a solution of the dKP hierarchy, i.e., satisfies 

d Tn C = {B n ,C}, (3.6) 
where a e denotes the principal symbol [7], the Poisson bracket is defined as 

= __-__, (37) 

and B n = (£)> now refers to powers of p. 

The dKP hierarchy can be also written in the zero-curvature form 

BE 8B 

^-^ + {B n ,B m } = 0. (3.8) 

In [7], from (2.3) and (2.4) (with L, B n , M and <9 n replaced by L e , B en , M e and ed Tn 
respectively ), it has been proved that ^(7) has the following WKB asymptotic expansion 
as e — > 

^) = ea;p[^S , (T,A) + 0(l)], e -> 0. (3.9) 

In order to take the quasiclassical limit of the constraint (3.2), we also need to find the 
asymptotic form of the adjoint Baker function. Similarly like the proof in [7], from (2.5), we 
can find has the following WKB asymptotic expansion 

ip*(^)=exp[~S(T,X)+0(l)], e^O. (3.10) 



From (2.3b) and (3.9), we obtain a hierarchy of conservation equations for the momentum 

as 
ax> 

dp _ dB n ( P ) 



function p = jrj?. 



dT n dX ' 

the compatibility of which, i.e., dT q T = 9T | T implies the dKP hierarchy (3.1 
Regarding 

ft (|) = A = ~ exp[ SiT ^ Xi) + a n + 0(e)], e - 0, (3.12a) 

^(7) = ^(7' A = ~ e ^[~ e ' A + ai2 + °( € )1' e ^ ' i = 1 »->^ ( 3 - 12b ) 
we will find that 

= e a - + ^[(ea)- 1 + (^g^ + 0( e ))(ea)- 2 + ((^g^) 2 + 0(e))( e a)- 3 + •••], e ^ 0. 
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Taking the principal symbol of both sides of (3.2), we have 

N 

C n = B n + Y,e a ^[p- 1 + S x (T,X z )p- 2 + S 2 x (T,X l )p~ 3 + ---] 

i=l 
N 

= B n + J2^- (3.13) 
where B n = (£™)>o and 

u . = e «ii+«« > p i = S x (T,X i ). (3.14) 

The constraint (3.13) is well known and is often called Zakharov reduction when n = 1 
[2, 9, 10, 12]. In [12], using the 3 method, the authors demonstrated that (3.13) is a result of 
symmetry constraint of the dKP hierarchy. Here we have shown that (3.13) can be obtained 
by the dispersionless limit of (2.6), i.e., the dispersionless limit of symmetry constraint for 
the Lax operator of the dispersionfull hierarchy leads to the symmetry constraint for the 
Lax function of the dispersionless hierarchy. From (3.3), (3.12), (3.14) and by a tedious 
computation, we obtain the following equations of hydrodynamical type 

Pi,T k = [B k (p)\p= Pi \x, (3.15a) 
d 

Vi,T k = [vi(-^B k (p))\ p=Pi }x, i = l,---,N. (3.15b) 

Taking the dispersionless limit of (3.4), we will get the constrained dKP hierarchy reduced 
by (3.13) 

(C n ) Tk = {B k ,C n }, (3.16a) 

Pi,r k = [B k (p)\p= Pi }x, (3.16b) 
d 

Vi, Tk = [vi(-7^B k (p))\ p=Pi \x, i = 1, -., N. (3.16c) 

Adding the term (i3fc)r n to the right hand side of (3.16a), or taking the dispersionless limit 
of (2.9) directly, we will obtain the dKP hierarchy with self-consistent sources as 

(B k ) Tn -(C n ) Tk + {B k ,C n } = 0, (3.17a) 

Pi,r k = [B k (p)\ P = Pi \x, (3.17b) 
d 

Vi, Tk = [vi(-^B k (p))\ p=Pi ]x, i = 1, -,N. (3.17c) 

If "T n " is viewed as the evolution variable, (3.16) may be regarded as the stationary system of 
(3.17) like the dispersionfull case. It is not difficult to prove that under (3.17b) and (3.17c), 
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the equation (3.17a) will be obtained by the compatibility of the following conservation 
equations 

p Tk = (B k (p)) x , (3.18a) 

N 

Pr n = (£ n (p))x = \B n {p) + ( 3 - 18b ) 



For example, when n = 3, k = 2, 



N 

Vi 



£ 3 = p 3 + 3U lP + 3U 2 + Y — B 2 =p 2 + 2U U 

~^P~Pi 

(3.17) becomes the dKP equation with self-consistent sources (U 2 is eliminated by U 2 ,x 
\U 1X and T = T 3 ,Y = T 2 ) 

N 



{U lT _ +J2 v i,x)x = jU 1>YY , (3.19a) 

i=l 



Pi ,Y = (p 2 i +2U 1 ) x , (3.19b) 

v ii Y = 2(v i p i )x, i = l,...,N. (3.19c) 

Under (3.19b) and (3.19c), the compatibility of the following conservation equations give 
rise to (3.19a) 

p Y = (p 2 + 2U 1 ) X = 2 PPx + 2U hX , (3.20a) 

N N N , x 

p T = (p 3 +3u lP +3u 2 +y2 —)x = 3p 2 px+3(u lP )x+3u 2 ,x+y2 — -E Vi{ ? x ~ p :?K 

j^p-pi j^p-pt [P-Pi) 2 

(3.20b) 

where U 2 ,x = 



4 Hodograph solutions 

Motivated by the dKP case, we would like to solve (3.17) by hodograph transformation 
provided the conservation equations (3.18). Following [5], one can consider the M-reductions 
of (3.18) so that the momentum function p and Vi, pi,i = 1, N depend only on a set of 
functions W =(Wi,...,W M ) with Wi = U\ and (W 1: ...,Wm) satisfy commuting flows 

dW dW 

Wr ^(W) w , n>2 (4.1) 
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where the M x M matrices A n are functions of (W±, Wm) only. In the following, we shall 
take the dKP equation with self-consistent sources (3.19) for example and show its solutions 
in the cases of M = 1 and M — 2. 
1. M = 1 

In this case, p = p(Ui), Vi = Vi(Ui), p, t = Pi(Ui) and (4.1) becomes 

U 1X = A(U 1 )U 1 , X , U 1:T = B(U 1 )U 1:X . (4.2) 

From (3.19b), (3.19c) and (3.20), we will get the following relations respectively 

dpi f A(Ui) . . 

dvi ,A(Ui) , dpi 

du^—2 — ti = Vi mi> (43b) 

dp ,A(U, 



dU x 



1J --P) = 1, (4.3c) 



7 , , jtt N dvi N ( dp dpi \ 

dp d(tt\ q 2 dp dp dU 2 sr^ Wl ^Wl ~ dul) (A oa\ 



i=l ^ ■ r ' i=l 

which implies 

5 = 3^ 1 + -A 2 -> — -. A = 2 — -. (4.4) 
1 A ^ dUi dUi y ! 

It is easy to verify that with (4.3) and (4.4), (4.2) are compatible. Making the hodograph 
transformations with the change of variables (X,Y,T) — > (Ui,Y,T) and X = X(Ui,Y,T), 
we will get the following hodograph equations for X, 

dX A dX „ 3 <9 (ifj 

i=l 

which can be easily integrated as follows 

o n , 

X + A(^)y + (3^ + -A^) 2 - £ ^-)T = FM, (4.6) 

i=i 1 

where F(Ui) is an arbitrary function of U\. 

If we choose A(Ui) = C = const, F(Ui) = (3Ui, from (4.6), (4.3a) and (4.3b), we get an 
implicit solution as 

Q N C 2 

X + C Y+ [3U, + -Cl - - + <U))-*]T = 0U U (4.7a) 

i=i 
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v i = c i [^f-2(U 1 + d i )]-K (4.7b) 



Pl = Y ± Vx _2( ^ 1 + di) ' i = 1 >-> JV > ( 47c ) 

where = 1,...,N are constants. If aq = d 2 — • • • — d>N an d Z)i=i c « = 0> (4.7a) 

degenerates to the solution of dKP equation [5] . 

If we assume Vi — v 2 — • • ■ — v N and J2?=iju^ = 3Z7i, i.e., = jjUi, and taking 
F{U\) = 0, by a direct computation, we will obtain an explicit solution of the dKP equation 
with self-consistent sources (3.19) as follows 

_ AY 2 - QTX ± AY VY 2 - 3TX 
1 ~ 225T2 ' 

3 o 
1 2N 1: 



Pi 



2y/2Ul, i = l,...,N. 




2. M = 2 

In this case we denote W 1 = U 1 , W 2 = W, then v t = Vi(Ui,W), pi = Pi(Ui,W) and 
p = p(Ui, W) with the commuting flows 

where A = (A^) and B = are 2x2 matrix functions of Lq and W. Requiring U ljX and W x 
are independent, (3.196), (3.19c) and (3.20) give rise to the following relations respectively 

dpi dpi dpi dpi 

dv^ dv^ _ d(v iPi ) d(v tPt ) 

*)A = ^# *) + (2,0), (4.9c) 



( dp dp_\n _ q„2/ dp dp \ , n ( d(Uip) dJUip) \ , o/gl/ 2 g£/ 2 \ 

+ V W 1 f ^» _ _dp_ _ dp^ dp dpi \ 

which implies A(Lq,VF) and W) must satisfy 



(4.9d) 



i=i 1 \Ai2 Z^i=i aw Ai 2 Z^i=i aw. 



N dvi 
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with A u = 2fg and A 12 = 2§£. 

For simplicity we assume = 0, i — 1, N. Then 

B = \a> + w -t,m: I= + ^ - i rfrf - 4 - 1 E («!) 

i=i j=i 

where the formula A 2 = (trA)A — (detA)I is used. 

With (4.11), the compatibility for (4.8) requires A to satisfy 

1 ddetA 



1 dtrA 



4 3 ^i=l at/i ' _ 

dUl / V 4 dUi 



Using the hodograph transformation changing the independent variables (X, Y, T) to (Lq, W, T) 
with X = X(U U W, T) and Y = Y{U U W, T), we get 



where §|^^y = X W Y T — X T Y W . It is not difficult to see that (4.13) has solutions in the form 

1 1 N d 

X + 3(C/i - -detA g^)T = F(U U W), (4.14a) 

Y + -(trA)T = G(U!,W), (4.14b) 

where we have required that Y u± and IV are independent, while F and G satisfy the linear 
equations 

"M = A ( Gw ) . (4.15) 



An example of solution is given by 

4 W ' 



^= , „J > ( 4 - 16 ) 



and 

Vi = CiUi, Pi = Y> i = 1 '"'' iV ' ( 4 - 17 ) 
with Cj,i = 1, • • • , TV are constants. 
(4.14) now becomes 

1 1 N 

X + 3[U 1 ~- A {W 2 - 4[q) C *\ T = F (Ui, W), (4.18a) 

i=i 
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Y + -2WT=G{U U W). (4.18b) 
From (4.15) and F WUl = F UlW) G must satisfy 

2Gu 1 + UiGjj 1 u 1 — AGww — 0. 

Taking G = -W and from (4.15), F = \W 2 - AUi and we obtain a solution of (3.19) as 
follows 

tj Y 2 | (El^T-X 

1_ 2(3T + 2)2 + 2(3T + 2) ' 1 ' 

Y 2 ( S T N c )T - X 

Vt = CiUl = Ct [" -+ tJ - - 1, (4.19b) 

[ 2(3T + 2) 2 2(3T + 2) J ' v ; 

W Y 

Pi = — = . i = l,...,N. (4.19c) 

F 2 3T + 2 v ; 

(4.19) is a global solution for T > — | and (4.19a) degenerates to the solution of the dKP 

equation [5] when Y^=i c « = 0- 
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